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Abstract 

I review a simple method, recently introduced to convert rheological 
compliance measurements into frequency-dependent moduli. New exper- 
imental data are presented, and the scientific implications of various data 
conversion methods discussed. 

This is an article about the sanctity of experimental data. 

As the readers of this Bulletin will know better than I, there are many dif- 
ferent ways to do rheology — many different experimental techniques by which 
to measure the mechanical properties of viscoelastic materials. Some tech- 
niques measure a compliance, others measure a modulus. Some measurements 
are most conveniently expressed as a function of time, others as a function 
of frequency. For instance, in a passive microrheology experiment, the time- 
dependent mean-square displacement of probe particles is measured as the 
particles perform Brownian motion within the viscoelastic medium. That in- 
creasing function of time is a measure of the fluid's time-dependent compliance 
[2 [21(3]. There is often no alternative to such microrheological measurements, 
as sample volumes may be very limited, due to expense, or difficulty in produc- 
ing or refining the samples [4] . The complance is also provided by a macroscopic 
creep experiment where a stress step is applied to a material, and its resulting 
strain is followed as time proceeds. Even where sufficient material exists to use 
a conventional rheometer, a creep experiment is often the only feasible probe of 
the slowest relaxation modes of high- viscosity materials such as drilling muds 
or soft glassy materials including pastes and slurries. The raw data generated 
by those experiments look quite different from the frequency-dependent storage 
and loss moduli (G'(u)) and G" (io)) measured in oscillatory rheometry. Nev- 
ertheless, so long as the measurements were performed in the linear regime, 
the two experiments in fact probe exactly the same physical properties of the 
material, so that the time-dependent compliance and the frequency-dependent 
modulus data contain the same information. 

This can be easily understood if one imagines performing oscillatory mea- 
surements at two different frequencies simulateously — a technique known as 
multi-wave rheology. Applying a high-frequency oscillatory stress superim- 
posed on top of a low-freuqency oscillatory stress causes the material to respond 
at both frequencies simultaneously, and the high- and low-frequency parts of its 
strain can be separately measured, giving its moduli at both frequencies. The 
experimenter could add in more frequencies, to measure more points on the 
spectrum at once. Every schoolboy who has watched an oscilloscope screen as 
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sine waves from different oscillators were added together knows that, by adding 
up the appropriate amounts of all the odd harmonic frequencies, one can gen- 
erate a square wave. A single stress step is just a square wave of infinitely low 
frequency, and is therefore equivalent to performing oscillatory measurements 
at all frequencies at once. 

Since the data from creep and oscillatory rheometry contain equivalent in- 
formation, it is equally acceptable to express their results by plotting a graph 
of compliance versus time, or of storage and loss moduli versus frequency. It is 
useful to express the results from both experiments in the same format, so that 
they can be directly compared. In practise, the frequency-dependent moduli 
have become the more standard rheological measure, so the need exists to con- 
vert creep compliance data into G'(u>) and G"(to), and therein lies a problem. 
How do we convert the data accurately, i.e. without introducing artefacts? 

In principle, the conversion is absolutely straightforward, since the frequency- 
dependent moduli are inversely proportional to the frequency-dependent com- 
pliance. We have measurements of the tame-dependent compliance, so we 
just have to identify how much of each frequency is contained in that time- 
dependent signal, i.e. Fourier-analyse it. There are plenty of standard tools 
available for Fourier- analysing a signal, but the compliance signal from a step- 
stress experiment has several features that over-tax the approximations inher- 
ent in those tools, liberally decorating the output with artefacts. We need to go 
back to the drawing board, and analyse closely what is involved in processing 
a signal, making sure not to violate basic principles of scientific method. 

The main problem is one inherent in many experiments. We want to process 
(in this case, to Fourier-analyse) a signal that is a function of time. Unfortu- 
nately, we do not have access to the true function of time since, in practice, the 
data have three shortcomings: the measurements begin some short but non- 
zero time after the stress step was applied (due to the finite time resolution 
of the rheometer); the measurements are discrete, not continuous (also due to 
the temporal resolution); and the data terminate after some finite duration. 
The mathematical recipe for a Fourier transform requires as input a function 
that is defined for all values of the time t. So we have a dilemma: what con- 
tinuous function of t should we substitute into the mathematics instead of the 
experimental data? The knee-jerk reaction to this question is one of horrified 
indignation: don't substitute an artificial function; just stick the experimental 
data into the Fourier-transformation formula. Sadly, the Fourier-transform for- 
mula demands a continuous function, so blindly substituting-in discrete data 
is equivalent to replacing a data-set such as that in Fig. [T^ by the function 
in Fig. [lb. Clearly, some interpolation and extrapolation is unavoidable, and 
we need to find a fitting function that will be truest to the data. Those who 
try to avoid choosing a particular function by simply substituting their data 
directly into the Fourier transform are unwittingly choosing to interpolate with 
a function such as that in Fig. [TJd. There could hardly be a worse choice. It 
would be similarly unwise to evade the choice by entering the data into some 
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Figure 1: (a) Sketch of an example data set from experimental measurements 
of a time-dependent compliance J. (b) Naively substituting the data directly 
into a Fourier integral (replacing the integral by a summation over the data 
points, J2 P Jp exp(iu> t p )) is tantamount to replacing the data by the function 
shown (a series of delta functions), and thus assuming that the compliance has 
a value of zero at all times except for those when a measurement was taken. 
That is obviously a bad assumption. 

proprietary spectrum analyser. The approximation scheme would then be hid- 
den from the scientist, but the responsibility for the resulting artefacts would 
remain his. Let us review the scientific validity of some other popular choices 
of fitting function for the compliance data. 

A common approach is to fit the data to a theoretical model of viscoelastic 
fluids, such as the generalized Maxwell model, in which a number of Maxwell 
modes are added together to model relaxations in the fluid that occur at a range 
of different rates [5]- Each mode is characterized by an amplitude (how much 
of that mode is present) and a time scale. Thus the experimental compliance 
data are used to establish the values of a set of fitting parameters — two 
per Maxwell mode — and the frequency-dependent moduli of that best-fit 
model can then be plotted. The philosophy of such a scheme is difficult to 
justify, since it blurs the distinction between theory and data. Should the 
resulting graphs of G'(uj) and G"(u>) be regarded as the results of a model 
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or an experiment? They can be legitimately interpreted as the output of a 
model, with the understanding that they therefore do not carry the weight of 
a true empirical observation. It would be unscientific to treat the results as 
measurements, since scientific methodology demands that experimental design 
must not allow data to be influenced by theoretical preconceptions. It must 
not be forgotten that experimental data, not theoretical models, are sacred. 

In the same way that sine waves of different frequencies can be added to- 
gether to reproduce any function (a Fourier transform), the same is true of 
decaying exponentials (a Laplace transform [5]); if one adds together suffi- 
ciently many decaying exponentials of appropriate amplitudes and rates, one 
can fit them to the experimental data with arbitrarily good accuracy. So a 
large series of exponentials need not be regarded as "Maxwell modes" derived 
from a theoretical model, but simply as a general smooth function with suf- 
ficiently many fitting parameters to give a good approximation of the data. 
Using that fitting function has the advantage that it is mathematically very 
easy to Fourier-transform, and thereby convert into graphs of G'(u>) and G"(u)). 
Philosophically, we are now of firmer ground. So long as the fitting function 
has enough flexibility (enough fitting parameters) to accurately reproduce the 
data set, it can be legitimately regarded as a true, unbiassed representation 
of the experimental data, so that the resulting graphs of G'(u>) and G"(oj) are 
scientifically valid experimental results. But how many fitting parameters are 
enough? How close a fit to the data is close enough? 

Intuition might suggests the following critereon. Data from the compliance 
experiment inevitably have some scatter, such as depicted in Fig. [5^. If a fitting 
function is so inflexible (like the dashed curve) that it can't be made to pass 
through the middle of that scatter, then it is clearly unsuitable, as it is a biased 
(inaccurate) representation of the data. If, on the other hand, by fitting many 
parameters, one can find a smooth fitting function that passes through the 
middle of the noise (like the solid curve in Fig. \2&) , then it is usually assumed 
to be good enough. I disagree. 

It is our sworn duty as scientists not to doctor data, but to report them 
accurately, and every graduate knows that the best way to treat experimental 
noise and ouliers is not to conceal them, but to present them for the reader's 
scruitiny. Therefore, I want to argue that the only acceptable fitting function 
is one that passes through every data point. Only then does the interpolating 
function truthfully report the results of the experiment; it is a true represen- 
tation of the sacred experimental data. It is important to keep in mind the 
purpose of this exercise in interpolation and extrapolation, which is to allow 
us to Fourier-analyse the data, not to smooth them. On our quest to find a 
suitable interpolation scheme, we must not be beguiled by shiny things like the 
lovely smooth fitting function in Fig. [2^. The assumption that the compliance 
should be smooth is a prejudice based on our preconceptions about fluids, but 
the experimental apparatus did not report a smooth function. We must put 
our preconceptions aside, and use the data. 
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Figure 2: (a) Sketch of an example data set of a time-dependent compliance 
measurements with noise, and two alternative smooth fitting functions, (b) A 
picewise-linear interpolation between the data point, and also piece- wise linear 
extrapolations. (Inset: An alternative interpolation scheme in which the curve 
is unnecessarily constrained to have continuous gradient at all points.) 

So, to turn the discrete set of data points into a continuous function, suit- 
able for Fourier analysis, we must join the dots, as in Fig. [5b. This still leaves 
us with an infinite choice of functions with which to connect the data points, 
and one might still be tempted to construct a locally smooth function by using, 
for instance, a piecewise cubic or a sine series. However, such functions have 
a tendency to overshoot, as depicted in the inset to Fig. [2Jd. By insisting un- 
necessarily that the gradients match where the curves connect, in some cases 
one forces the curves into contrived shapes that actually lie far from the data 
set. As well as passing through every data point, and thereby telling the whole 
truth, the chosen fitting function must tell nothing but the truth, i.e. it must 
have no unwarranted behaviour in between the data points, as that would in- 
troduce spurious new information that was not present in the data. The only 
function that reliably introduces no new features between two data points is a 
straight line segment. So the most truthful interpolation of the data is given 
by a piecewise-linear function, as sketched in Fig. [5b. 

To turn a finite data set into a function defined for all t, one must not 
only interpolate but also extrapolate to t — and t = oo, which we do with 
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two additional straight-line segments. If the creep experiment was performed 
for sufficiently long to reach the terminal regime, where all elastic memory is 
lost, and the fluid flows with its long-time viscosity 77, then the straight-line 
extrapolation to t = 00 is uncontroversial. The final line-segment is given a 
gradient I/77, found by esimating the asymptote to the data by eye. Note 
that the human eye is the only good way to judge an asymptote to noisy or 
unpredictable data, as automated schemes are easily fooled. Similarly, the 
value Jo of the compliance at t — must be carefully judged from a linear or 
log-linear plot. If, on the other hand, the data do not reach the terminal regime, 
then no amount of analysis can recover the frequency-dependent moduli from 
the incomplete experiment. 

It turns out to be straightforward to write down the exact Fourier transform 
of the piecewisc-linear fit, yielding a simple formula. Plugging the data points 
into the formula is child's play, so that some standard laborious methods for ap- 
proximate conversion are rendered redundant [TJ [7J E] . The simple formula can 
be found in Ref. [8j, or can be downloaded for various standard data process- 
ing packages from www. pcf . leeds . ac . uk/research/h.ighlight/view/4 but 
is also reproduced here 



in terms of the two parameters Jo and T], and the N data points (if*, Jk)- This 
formula will return an answer for any value of frequency lo that is substituted 
into it. But, of course, the data only contain information about a limited range 
of frequencies. Specifically, the frequency window depends on the resolution (of 
order t-y, the time of the first data point) and duration £jv of the data set, so that 
the formula returns real, scientifically valid information for frequencies in the 
range i^ 1 < u < t± , and artefact-dominated rubbish for frequencies outside 
that range. To be clear, outside that frequency range, the function is strongly 
influenced by the artificial features of the fit: the straight-line extrapolation 
to infinity and the sharp corners of the jagged function, whereas all features 
(including noise) that appear on the modulus graphs within the operational 
frequency range are genuine features of the data. 

To demonstrate the efficacy of the conversion formula, new data are pre- 
sented here (courtesy of Chirag Kalelkar, Complex Fluids and Polymer Engi- 
neering Group, National Chemical Laboratory, Pune, India) for a PAS melt. 
Figure [3] shows the compliance measured using CP-25 geometry on an Anton 
Paar MCR301 rheometer, at a stress of 200 Pa, averaged over 10 independent 
creep compliance runs with a duration of 600 s, combined with three measure- 
ments of duration 1 hour. It is clear from the figure that the duration was suffi- 
cient to reach the terminal regime, where the complance asymptotes to a linear 
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Figure 3: Creep compliance J versus time t for a polyalkylalkenylsiloxane 
(PAS) melt (molecular weight M w = 5 x 10 5 g/mol, polydispersity index 
M w /M n = 1.83) measured at room temperature using 25 mm cone-plate geom- 
etry on a stress-controlled Anton Paar MCR301 rheometer, at a stress of 200 Pa, 
averaged over 10 independent creep compliance runs up to time t = 600 s. 
Strain rate measurements averaged from three longer runs, also at 200 Pa, 
at lower temporal resolution were used to extend the range to t = 3600 s. 
The linear extrapolation beyond 3600 s is also shown, using the parameter 
fj = 103600 Pa s. Inset: close-up on the early data points, shown on linear 
axes, to allow the zero-time extrapolation of J to be estimated. Data courtesy 
of Chirag Kalelkar, Complex Fluids and Polymer Engineering Group, National 
Chemical Laboratory, Pune, India 

function (i.e. the Newtonian regime), with viscosity r\ = 103600 ± 100 Pa s, as 
demonstrated by the linear extrapolation shown in Fig. [31 The inset to Fig. [3] 
shows the first few compliance data points on linear axes, allowing us to esti- 
mate the zero-time value as Jq — (1.12±0.09) x 10~ 5 Pa -1 . Figure 0] shows the 
results of plugging the compliance data into the formula for G'(u>) and G"(oj) 
(continuous curves). Also, superposed are data from oscillatory measurements 
(frequency sweep) on the same material. 

Note the striking agreement between the two measurements, confirming 
that the different experimental techniques probe the same material proper- 
ties. Note also that the converted creep measurement just reaches the terminal 
regime (confirming that the compliance measurements were just long enough) , 
where G oc uj 2 and G" oc w, at frequencies that are impractically low for oscil- 
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Figure 4: Storage and loss moduli as functions of angular frequency. The 
continuous curves were found by plugging the compliance data of Fig. [3] into 
the formula for G'(u>) and G"(u>). They are plotted only within the valid 
frequency range tjj 1 < u> < , in this case 2.8 x 10 _4 s _1 < u) < 20 s -1 . Also, 
shown for comparison are data (courtesy of Chirag Kalelkar) from oscillatory 
measurements (frequency sweep) on the same material (G' grey circles, G" grey 
squares). 



latory measurements. Since Fig. U displays only the relevant frequency range 
t^ 1 < lu < t^ 1 , this low- frequency behaviour is not an artefact of the straight- 
line extrapolation; it is a result dominated by the late-time data that have 
reached the terminal regime in the creep experiment, while the straight-line 
extrapolation controls features at lower frequencies than those plotted in the 
figure. Only if a wildy inaccurate value of the viscosity were used, that was at 
odds with the data, thus introducing an abrupt corner into the overall shape 
of the fitting function, would the resulting artefacts be so large as to leak into 
the operational frequency window. 

The new conversion formula has already been applied to a number of differ- 
ent experimental systems in addition to the PAS melt reported above, including 
a polyisoprene melt in conc-and-plate geometry [5] , magnetic microrheological 
measurements of F-actin solutions (a semi-felxible biopolymer) [9] , and semidi- 
lute polyacrylamide solutions measured by both magnetic microrheology [5] 
and a new protocol for obtaining microrheological creep measurements using 
optical tweezers |10j . The technique has even been applied to non-equilibrium 
samples of aqueous wormlike micellar solutions [llj . since it does not require 
the sample to be at equilibrium. So long as the sample's moduli evolve only 
slowly — on a much longer time-scale that the measurement — the conversion 
remains valid. This makes creep the measurement of choice for slowly evolving 
non-equilibrium samples, as it can be performed more quickly than the equiva- 
lent oscillatory frequency sweep. Of course, the formula should not be applied 
to non-linear rheology or to materials that are aging rapidly on the timescale 
of the measurements. In both of those cases the dynamic moduli are anyway 
not uniquely defined, as their definitions must then take account of coupling 
between harmonics [T2j and of absolute time [13]. As the conversion method 
was published only a matter of months ago [8] , the list of relevant literature is 
rather short, but growing rapidly due to the ease with which the formula can 
be applied. 
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